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ABSTRACT 

The densities in the outer regions of clusters of galaxies are very low, and the coUisional timescales 
are very long. As a result, heavy elements will be under- ionized after they have passed through the 
accretion shock. We have studied systematically the effects of non-equilibrium ionization for relaxed 
clusters in the ACDM cosmology using one-dimensional hydrodynamic simulations. We found that 
non-equilibrium ionization effects do not depend on cluster mass but depend strongly on rcdshift 
which can be understood by self-similar scaling arguments. The effects are stronger for clusters at 
lower redshifts. We present X-ray signatures such as surface brightness profiles and emission lines 
in detail for a massive cluster at low redshift. In general, soft emission (0.3-1.0 keV) is enhanced 
significantly by under-ionization, and the enhancement can be nearly an order of magnitude near 
the shock radius. The most prominent non-equilibrium ionization signature we found is the O VII 
and O VIII line ratio. The ratios for non-equilibrium ionization and coUisional ionization equilibrium 
models arc different by more than an order of magnitude at radii beyond half of the shock radius. 
These non-equilibrium ionization signatures are equally strong for models with different non-adiabatic 
shock electron heating efficiencies. We have also calculated the detectability of the O VII and O VIII 
lines with the future International X-ray Observatory (IXO). Depending on the line ratio measured, 
we conclude that an exposure of ~130-380 ksec on a moderate-redshift, massive regular cluster with 
the X-ray Microcalorimeter Spectrometer (XMS) on the IXO will be sufficient to provide a strong test 
for the non-equilibrium ionization model. 

Subject headings: galaxies: clusters: general — hydrodynamics — intergalactic medium — large-scale 
structure of universe — shock waves — X-rays: galaxies: clusters 



I 1. INTRODUCTION 

Clusters of galaxies are very sensitive probes 
• to cosmological par ameters (e.g., lAllen et "aD 120081 : 
. iVikhlinin et al.l [20091 ) . Systematic uncertainties in pre- 
cision cosmology using galaxy clusters can be minimized 
by restricting the sample of clusters to the highest de- 
gree of dynamical relaxation, and hence studying relaxed 
clusters is particularly important. In addition, the outer 
I envelopes of clusters have been thought to be less sub- 
jected to complicated physics such as active galactic nu- 
cleus feedback, and hence these outer regions may pro- 
vide better cosmological probes. 

Clusters are believed to be formed by the continual 
merging and accretion of material from the surrounding 
large-scale structure. Based on high resolution A^-body 
simulations, it has been found that cluster growth can 
be divided into an early fast accretion phase dominated 
by major mergers, and a late slow phase dominated by 
smooth accr etion of background materials and man y mi- 
nor mergers (jWechsler et al.|[200l[Zhao et al.|[2009l ). Ac- 
cretion shocks (or virial shocks) are unambiguous predic- 
tions of cosmological hydrodynamic simulations. For the 
most relaxed clusters with roughly spherical morphol- 
ogy in the outer regions, these simulations show that 
large amounts of material are accreted through filamcn- 
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tary structures in some particular directions, and more 
spherically symmetrically in other directions. 

Most of the simulations assume the intracluster 
medium (ICM) is in coUisional equilibrium. However, be- 
cause of the very low density in the cluster outer regions 
^20(0), the Coulomb coUisional and coUisional ion- 
ization timescales are comparable to the age of the clus- 
ter. It has been pointed out that electrons and ions there 
may be in non-equipartition (|Fox &: Lo cb"1997': 'Takizawa| 
[T999l[2000l:rRudd fc Nagaill2009l: IWong fc Sarazin.2009lL 
and also the ionization state may not be in col- 
lisional ionization equilibrium ("Yoshikawa e t al.l 120031 : 
ICen fc Fan j [20061: TYoshikawa & Sasaki 2006). If these 
non-equilibrium processes are not properly taken into ac- 
count, the measured properties may be biased in these 
regions. 

Studying the outer regions of relaxed clusters not only 
is valuable in understanding the accretion physics and 
to test the assumptions concerning plasma physics near 
the shock regions, but it is also very important to test 
structure formation theory and to constrain the system- 
atic uncertainties in clusters to be used as precision cos- 
mological probes. Before the launch of the Suzaku X- 
ray observatory, physical properties such as tempera- 
ture have never been constrained with confidence be- 
yond roughly one-half of the shock radius. Most of 

-Ra is the radius within which the mean total mass density 
of the cluster is A times the critical density of the universe. The 
virial radius iJvir is defined as a radius within which the cluster is 
virialized. For the Einstein-de Sitter universe, iJvir ~ ^il78i while 
for the standard ACDM Universe, iJvir ~ ^95- 
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our understanding of these regions is still based on the- 
oretical models and hydrodynamic simulations. Re- 
cently, observations by Suzaku have constrained tem- 
peratures up to about half of the shock radius for a 
few clusters for the firs t time to bett er than a fa ctor of 
~ 2 (IGeorge et al.ll2009l: IReiprich et al . 2009; Basu et al.l 
I201C : IHoshino et al.ll201Cll) . While the uncertainties are 
still large, there is evidence that the electron pressure in 
cluster outer regions may be lower than that predicted by 
numerical simulations assuming collisional equilibrium. 
Observations of secondary cosmic microwave background 
anisotropics with the South Pole Telescope (SPT) and 
the Wilkinson Microwave Anisotropy Probe ( WMAP) 
7 year data also support these results (|Komatsu et al.l 
120101 iLueker et aIll2010D . These observational signa- 
tures are consistent with electrons and ions in non- 
equipartition, although it is also possible that the hy- 
drodynamic simulations may simply overestimate the 
gas pressure. Another possibility is that heat conduc- 
tion outside the c luster may be reducing the gas pres- 
sure (|Loebll2002t) . Recently, iWong et all H2O1O0 have 
shown that cosmological parameters will be biased if non- 
equilibrium effects (in particular, non-equipartition) are 
not properly taken into account. 

More observations of the outer regions of clusters are 
being done or analyzed. In the future, the proposed In- 
ternational X-ray Observatory (/XO jj will have the sen- 
sitivity to constrain cluster properties out to the shock 
radius. Thus, a detailed study of the physics in the outer 
regions of clusters would be useful. In particular, IXO 
will have sufficient spectral resolution to resolve many 
important X-ray lines, and hence the ionization state of 
the p lasma can be det e rmine d. 

In iWong fc SarazinI (|2009h . we have studied in de- 
tail the X-ray signature of non-equipartition effects 
in cluster accretion shock regions. In this paper, 
we extend our study to include non-equilibrium ion- 
ization in our calculations. Non-equilibrium ioniza- 
tion calculations have been considered in a number of 
cosmological simulations to study the very low den- 
sity warm-hot inter galactic medium (WHIM) sur round- 
ing g a laxy clusters (lYoshikawa et 31.1 .2003: C en fc Fangl 
120061: lYoshikawa fc Sasakil \2QQ^ . At galaxy cluster 
scales, similar non-equilibrium i onization calculations 
have focuse d on merging clusters ()Akahori fc Yoshikawal 
120081 120101 ). These studies generally agree that in the 
low density ICM and the WHIM, there are significant 
deviations from ionization equilibrium, and that the ef- 
fects on the X-ray emission lines are strong. In this work, 
we focus on the X-ray signatures of non-equilibrium ion- 
ization in the accretion shock regions of relaxed clusters. 
We also discuss the detectability of non-equilibrium ion- 
ization effects with the future IXO. 

The paper is organized as follows. Section [2] describes 
in detail the physical models and techniques to calcu- 
late non-equilibrium ion fractions and X-ray observables, 
which includes the hydrodynamic models we used (Sec- 
tion [2T]), and the ionization and spectral calculational 
techniques (Sections l2.2H2.3p . The overall dependence 
the non-equilibrium ionization effects on the mass and 
redshift of the cluster is presented in Section [3] Calcu- 
lated non-equilibrium ionization signatures are presented 



in Section 21 which includes descriptions of particular 
models we used to pres ent the results C Section |4.1|) . X- 
ray spectra fSection l4.2[) . surface brightness profiles (Sec- 
tion 231), ^-nd the ratio of intensities of O VII and O VIII 
lines (Section 231). We discuss the detectability of O VII 
and O VIII lines and non-equilibrium ionization diagnos- 
tics with IXO in Section |5l Section IHl gives the discussion 
and conclusions. Throughout the paper, we assume a 
Hubble constant of Hq = 71.9 ft.71.9 km Mpc~^ with 
^171. 9 = 1, a total matter density parameter of ilM,o — 
0.258, a dark energy density parameter of Ha — 0.742, 
and a cluster gas fraction of /gas — ^b/^M = 0.17, where 
Oft is the baryon density parameter for our cluster models 
in the standard ACDM cosmologj0. The clusters have a 
hydrogen mass fraction X = 76% for the ICM. 

2. PHYSICAL MODELS AND CALCULATIONAL 
TECHNIQUES 

2.1. Hydrodynamic Models and Electron Temperature 
Structures 

We calculate the X-ray emission spectrum from 
the outer regions of cluster of galaxies using one- 
dime nsional hydrodynamic simulations we have devel- 
oped ()Wong fc Sarazinll2009D . Radiative cooling is negli- 
gible in cluster outer regions so that it will not affect the 
dynamics of the plasma. This non-radiative condition 
allows us to calculate the hydrodynamics and then the 
ionization structure and X-ray emission separately. The 
hydrodynamic models simulate the accretion of back- 
ground material from the surrounding regions onto clus- 
ters through accretion shocks in the ACDM cosmology. 
The calculated hydrodynamic variables (e.g., density and 
temperature profiles) in the cluster outer regions are con- 
sistent with those calculated by three-dimensional simu- 
lations, and a detailed disc ussion on our hydrodyn amic 
simulations can be found in IWong fc SarazinI ()2009[ ). 

The hydrodynamic simulations were done using Eule- 
rian coordinates. In calculating the time-dependent ion- 
ization for the hydrodynamic models, we need to follow 
the Lagrangian history of each fluid element. In one- 
dimensional spherical symmetric systems, it is convenient 
to transform from the Eulerian coordinates to the La- 
grangian coordinates by using the interior gas mass as 
the independent variable 



mg(r) — An 



'dr. 



(1) 



where r is the radius, and pg is the gas density. To do 
this, we first determined the values of mg{r) for each 
of the grid zones in th e hydrodynamical simulations in 
IWong fc SarazinI ()2009l ) at the final redshift of z = 0. 
The values of the density and temperature were also de- 
termined for this final time step for each gas element. 
Then, for each earlier time step, the values of gas den- 
sity and temperature within the shock radius were de- 
termined by linear interpolation between the nearest two 
values of the interior gas mass on the grid at that time. 
When the values of rng{r) fall between the discontinu- 
ous values of the shocked and preshocked elements, we 
assume the material to be preshock. This slightly un- 
derestimates the ionization timescale parameter defined 



http:/ /ixo. gsfc.nasa.gov/ 



^ |http: / /lambda.gsfc.nasa.gov/product / map / dr3/parameters-summary.cfm| 
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below, but this effect is small as the grid zones in the 
simulations are closely spaced. 

The coUisional ionization and recombination rates and 
the excitation of X-ray emission depend on the electron 
temperature (Tg). Because the accretion shock may pri- 
marily heat ions instead of electrons due to the mass 
difference and the long Coulomb coUisional timescale be- 
tween electrons and ions, the electron temperature may 
not be the same as the ion temperatu r e near the accre- 
tion shock regions ijFox fc Loebl Il997t iWong fc SarazinI 
Il009). The de gree of non-equipartition depends on the 
non-adiabatic shock electron heating ef ficiency, which is 
defined as /? = ATe^non-ad/ATnon-ad in iWong fc SarazinI 
(|2009f ). where ATg^non-ad and ATnon-ad are the changes in 
electron temperature and average thermodynamic tem- 
perature due to non-adiabatic heating at the shock, re- 
spectively. While there are no observational constraints 
on P in galaxy cluster accretion shocks, observations in 
supernova remnants suggest that /? <C 1 for shocks with 
Mach number s similar to those in galaxy cluster ac c retion 
shocks (Ghav amian et al.ll2007[) . IWong fc SarazinI (|2009f ) 
have calculated electron temperature profiles for models 
with a very low electron heating efficiency (/? = 1/1800), 
an intermediate electron heating efficiency (/3 = 0.5), and 
an equipartition model (/3 = 1). The electron tempera- 
ture profiles of these models with different values of /? 
are used to calculate the ionization fractions and X-ray 
emission in this paper. 

The ionization timescale parameter of each fiuid ele- 
ment is defined as 

/to 
Ue dt , (2) 

where rig is the electron number density, tg is the time 
when the fluid element was shocked, and to is the time at 
the observed redshift. The ionization timescale param- 
eters for cluster models with accreted masses of Mgh = 
0.77,1.53, and 3.06 x lO^^ Mq are shown in Figure [B 
Most ions of astrophysical interests will not achieve ion- 
ization equilibrium for the io nization timescale parame- 
ters < 10^2 (,m-3 s (iSmith fc Hughes 201(1). 

2.2. Ionization Calculations 

In order to calculate the X-ray emission spectrum for 
non-equilibrium ionization plasma, the ionization state 
of each fiuid element has to be calculated. We assume H 
and He are fully ionized. 

In this paper, we consider only coUisional ionization 
processes, and ignore photoionization. To justify this, 
we compare the photoionization and coUisional ioniza- 
tion rate for the ions we are interested in. The photoion- 
ization rate of an ion i is given by 

Jo 

where fTi(z/) is the photoionization cross section of ion i at 
frequency v, and F{i') is the ionizing fiux. Near the clus- 
ter accretion shock, the dominant photoionization source 
is the UV background. We approximate the UV ionizing 
flux at z = as a power law F{v) = Fq{v/vo)^^ with 
fn = 10^2 e rg cm'^ s'^ Hz, vq = lO^^-^ Hz, and T = 1.3 
(|Haardt fc Madau,1996i) . For example, for the O VII and 



O VIII ions we are most interested in throughout this pa- 
per, adopting the p hotoionization cross sections given by 
iVerner et all (|l996l ) gives i?Pi(0 VII) = 5.2 x lO-^^ s'^ 
and i?Pi(0 VIII) = 1.7 x lO'^^ s^^. The slowest col- 
lisional ionization rate of the oxygen ions is of the or- 
der of 10^^^(ne/10~^ cm~'^) s~^ for typical temperatures 
> 1 keV) and densitie s in the outer regions of clusters 
Smith fc Hughes! 120101 ) . This is nearly two orders of 
magnitude higher than the photoionization rates of O VII 
and O VIII ions. Even at z = 3 w hen the UV back- 
groun d is roughly 80 times stronger (|Haardt fc Madaul 
119961 ). the densities in the outer regions of clusters will 
be higher by a factor of ^ (1 + = 64, and hence 
coUisional ionization rates will increase by a similar fac- 
tor. CoUisional ionization still dominates over pho- 
toionization. For heavier elements, the photoionization 
rates will be even smaller, and the coUisional ioniza- 
tion rates for ions up to those of Ni are all higher tha. n 
- 10-i^(ne/10-5 cm-3) (ISmith fc HughesI IMOl . 
Therefore, we assume that photoionization is not impor- 
tant in our calculations. 

For ionization dominated by coUisional processes at low 
densities, the ionization states for each of the ions of a 
given element X are governed by 

^ ^ne{Ci-l{Te)fi-l +ai{Te)fi+l 

-[a(re)-Ha.-i(re)]/j, (4) 

where fi = n{X^'^^^)/n{X) is the ionization fraction of 
the ion i with charge +i — l, n(X+*~^) is the number den- 
sity of that ion, n{X) is the total number density of the 
clement X, and Ci{Te) and ai(Te) are the coefficients of 
coUisional ionization out of and recombination into the 
ion i, respectively. In solving Equation we use an 
eigenfunction technique which is based on the algorithm 
developed bv lHughes fc Helfandl (jl985D , and the method 
is desc ribed in detail in Appendix A of iBorkowski et aD 
(| 1994( 1. The eleven heavy elements C, N, O, Ne, Mg, 
Si, S, Ar, Ca, Fe, and Ni are included in our calcula- 
tions. The eigenvalues and eigenvectors, which are re- 
lated to the ionization and recombination rates, used to 
solve Equation (|4|) are taken from the latest version of 
the nei version 2 model in XSPECQ ( version 12.6.0). The 
atomic physics used to calculated the eigenvalues and 
eigenv ectors in XSPEC are discussed in IBorkowski et al.l 
(|2001l ). and the ionization fractions in the latest version 
of XSPEC are calculated using the updat e d diel ectronic 
recombination rates from iMazzotta et al.l (|1998l ) . All of 
the eleven heavy elements are assumed to be neutral ini- 
tially when solving Equation which is a good approx- 
imation as long as the ionization states of the preshock 
plasma are much lower than that of the postshock p lasma 
(e.g.. IBorkowski et al.l[T99l [200TI : IjTSilll IMOffl . We 
have also tested this by assuming all the ionization states 
in the preshock regions are in coUisional ionization equi- 
librium initially at different temperatures of 10^, 10^ and 
10^ K. We confirmed that the final ionization fractions 
for O VII and O VIII are essentially the same. 

2.3. X-ray Emission Calculations 
^ http:/ /heasarc. nasa.gov/xanadu/xspec/ 
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Once the ionization states are calculated by solving 
Equation (jj]) , the X-ray emission spectrum could be cal- 
culated by using available plasma emission codes such 
as the Raymond- Smith code (R aymond & Smith 1977), 
the SPE X code (IKaastra et ali Il996i). and the APEC 
code (.Smith et al .| |2001f ). ^ chose to use a version 
of the APEC code as implemented in the nei version 
2 model in XSPEC to calculate the X-ray emission spec- 
trum. The routine uses the Astrophysical Plasma Emis- 
sion Database^ (APED) to calculate the resulting spec- 
trum. Thus, the atomic data we used to calculate the 
ionization fractions and spectrum are mutually consis- 
tent, both coming from the same nei model version 2.0 
in XSPEC. The pubhcly available line list APEC.nei.vll 
was used throughout the paper. The inner-shell pro- 
cesses are missing in this nei version 2. We have com- 
pared our results to an updated line list which includes 
the inner-shell processes but is not yet publicly released 
(K. Borkowski, private communications), and we found 
that the difference is less than 10% which is smaller than 
the ^ 30% uncertainties in the atomic physics of the X- 
ray plasma code. 

We assume the abundance to be a typical value for 
galaxy clusters, which is 0.3 of solar for all models, and 
use th e solar abundance tables of lAnders fc Grevess^ 
([198^). 

We calculated the X-ray emissivity in photons per unit 
time per unit volume per unit energy, e^;, for each fluid 
element in our hydrodynamic models at redshift zero, 
which can be expressed in terms of an emissivity func- 
tion, Aij, by HSarazinlll986l) 



(-E = As n^Up , 



(5) 



where is the proton density. The emissivity function 
A^; depends on the ionization fractions and the electron 
temperature, but is independent of the gas density. The 
projected spectrum in the rest frame is given by integrat- 
ing t he X-ray emissivity along the line of sight (|Sarazinl 
[1981 



(6) 



where / is the distance along the line of sight. The broad- 
band rest frame surface brightness in energy per unit 
time per unit area is then given by 



IsEdE, 



(7) 



where the integral is across the energy band of interest. 

3. DEPENDENCE OF NON-EQUILIBRIUM IONIZATION ON 
CLUSTER MASS AND REDSHIFT 

We consider how the degree on non-equilibrium ioniza- 
tion of a cluster depends on its mass and redshift. We 
first consider how the cluster parameters which affect the 
ionization depend on mass and redshift. In order to com- 
pare equivalent locations in the different clusters, we es- 
timate the ionization parameters at a given over-density 
radius Ra- We consider a simple self-similar scaling ar- 
gument for galaxy clusters. 

® http:/ /cxc. harvard.edu/atomdb/ 



By definition, the average total density within the ra- 
dius R/\ is 



PA = Apcrit(2) 



3A 
SttG 



(8) 



where Pciit{z) is the critical density at redshift z, and 
H{z) is the Hubble constant at redshift z. The quantity 
E{z) = H{z)/Ho, so that 

E^(z) = [17m(1 + + nnil + zf + f^A 

+ (1 - -^R- ^a) (1 + zf] 



^^Mil + zf + \-nM. (9) 

The final expression in Equation ^ follows from the 
assumption that the universe is flat [VLm + + f^A — 1) 
and the fact that the radiation density parameter is small 
(rj/j ^ 1) in the present-day universe. The gas density 
at Ra is roughly 

Pgas (i?A) « /gasPA (X A /gasS'(z) . (10) 

The inflow timescale is roughly 



<«(GpA)-'^'ociff(0)A-i/2[ij(2)]- 



(11) 



where i/f(0) = 1/Hq is the Hubble time at z = 0. Thus, 
the ionization timescale parameter t varies as 

T cx Ai/2 f^^,E{z) . (12) 

Thus, the value of r at a fixed characteristic radius should 
be nearly independent of the cluster mass but will in- 
crease significantly with redshift. Figure [T] shows that 
the variation of t with radius is fairly self-similar for 
clusters with differing masses. When these curves are 
scaled to a fixed cluster characteristic radius, they are 
very nearly identical. 

The ionization state of the gas depends on r and the 
collisional ionization and recombination rates. For an 
under-ionized plasma, the collisional ionization rates are 
more important, so that the ionization state should de- 
pend mainly on Ci(Te)r. For an under-ionized plasma 
where the electron temperature is greater than the ion- 
ization potential of the relevant ions, Ci(Te) varies slowly, 

1 /2 

Ci(Te) (X Te ■ We have confirmed that this dependence 
fits the temperature dependence of the ionization rate of 
O VII over the interesting temperature range {kT^ — 1- 
5 keV). Below, we show that the ratio of O VIII to O VII 
lines is the best diagnostic of departures from ionization 
equilibrium (Section 14. 4|) . 

We assume that the electron temperature increases in 
proportion to the mean temperature Ta at the radius 
Ra- The mass within Ra is Ma = (47r/3)i?^pA, so that 
the radius Ra is given by 



Ra = 



2G 



1/3 



-2/3 



(13) 



The condition of hydrostatic equilibrium or the shock 
jump condition at the accretion shock imply that kTA ~ 
GMa/ Ra- Thus, the gas temperature varies as 



Ta oc il4/'[£;(z)]2/3. 

This implies that 

a{T,) o, Ml^'[Eiz)]'/\ 



(14) 



(15) 
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Combining Equations and (fT51) gives 

C,{T,)t M]l\E{z)f\ (16) 

Equations \12\ and (fT6|) suggest that the ionization state 
should depend only very weakly on the cluster mass, 
but should depend strongly on cluster redshift. The in- 
crease of r and Ci(Te)r with z implies that collisional 
ionization will be faster at high redshifts, and hence non- 
equilibrium ionization will be most important in low red- 
shift clusters. 

Later, we will show that the ratio of O VIII to O VII 
lines is the best diagnostic of departures from ioniza- 
tion equilibrium (Section 14. 4p . Figure [2] shows the di- 
mensionless ionization parameter Ci{Te)T for O VII as 
a function of the scaled radius (r/i?sh) for clusters with 
different masses at z = 0. Here, i?sh is the radius of 
the cluster accretion shock. Models with /3 = 1/1800 
(non-equipartition) and /3 = 1 (equipartition) are shown 
in thick and thin lines, respectively. All the curves for 
the same /3 but with different masses nearly overlap each 
others. This confirms that the non-equilibrium ionization 
effect at the same characteristic radius is nearly indepen- 
dent of mass. 

Figure [3] shows the dimensionless ionization parameter 
as a function of the scaled radius for cluster models at dif- 
ferent redshifts. The cluster model with a total accreted 
mass Msh = 1.53 x M© at z = is used. Models 
with (3 = 1/1800 (non-equipartition) and /3 = 1 (equipar- 
tition) are shown in thick and thin lines, respectively. In 
contrast to the dependence in mass, we find that there is 
significant evolution with redshift. The non-equilibrium 
ionization effect should be strongest for low-redshift clus- 
ters. 

The preceding arguments assumed that departures 
from ionization equilibrium could be assessed through 
the variation of Ci{Te)T. To test this explicitly and 
to determine if the preceding arguments apply to spec- 
tral diagnostics for non-equilibrium ionization, we study 
the effects of non-equilibrium ionization on O VII and 
O VIII ion fractions using numerical simulations. We 
calculate the non-equilibrium ionization bias, b, of the 
ionization fraction ratios /(O VIII)//(0 VII) for the non- 
equilibrium ionization (NEI) and collisional ionization 
equilibrium (CIE) models with the same electron heat- 
ing efficiency /3. The non-equilibrium ionization bias is 
defined as 

_ [/(o viii)//(o , . 

- [/(O VIII)//(0 VII)]ciE ■ ^ ^ 

Figure |4] shows the non-equilibrium ionization bias as a 
function of the scaled radius for clusters with different 
masses. Models with /3 = 1/1800 (non-equipartition) 
and /3 = 1 (equipartition) are shown in thick and thin 
lines, respectively. The nearly self-similar curves justify 
the semi-analytical argument that non-equilibrium ion- 
ization effect at the same characteristic radius is nearly 
independent of mass. Note that the non-equilibrium ion- 
ization effect on the /(O VIII)//(0 VII) ratio is only 
significant for the outer 10% of the shock radius. How- 
ever, the projected O VII and O VIII line emission will 
be significantly affected even at a radius as small as one- 
fourth of the shock radius. This is because the projected 
emission at the inner radius can be dominated by the line 



emission from the under-ionized outer shell (Section 14.41 
below) . 

Figure [5] shows the non-equilibrium ionization bias as 
a function of the scaled radius for cluster models at dif- 
ferent redshifts. The cluster model with a total accreted 
mass = 1.53 X 10^^ Mq at z = is used. Models with 
13 = 1/1800 (non-equipartition) and (3 — 1 (equiparti- 
tion) are shown in thick and thin lines, respectively. The 
non-equilibrium ionization effect is larger for clusters at 
lower redshifts which agrees with the semi-analytical ar- 
gument given above. Note that at a redshift higher than 
1, only a thin shell with a width of less than 5% of the 
shock radius is in non-equilibrium ionization compared 
to the wider shell with a width of ^ 10% of the shock 
radius for zero redshift clusters. 

4. NON-EQUILIBRIUM IONIZATION SIGNATURES 
4.1. Models Used to Calculate Spectra 

Massive clusters at low redshifts are ideal candidates to 
study the non-equilibrium ionization effects in the outer 
regions since the departures from ionization equilibrium 
are larger at low redshift and nearly independent of clus- 
ter mass (Section [3]). Clusters with higher masses are 
more luminous in X-rays, and hence the spectral signa- 
tures are easier to detect. In the following, we present 
X-ray spectra for the hydrodynamic cluster model with 
an accr eted mass of Msh = 1.53 x 10^^ Mq at z = calcu- 
lated in lWong fc Saraziiil ()2009n . This model represents a 
typical massive cluster in the present-day universe. The 
shock radius is i?sh — 4.22 Mpc for this model. The 
virial radius is i?vir — R95 — 2.75 Mpc, and the total 
mass within Rg^ is M95 — 1.19 x IO^^Mq. Another com- 
monly used radius and mass are i?2oo — 1-99 Mpc and 
M200 = 9.50 X 10"Mq. 

We calculate spectra for three different values of the 
shock electron heating efficiency /3 = 1/1800, 0.5, and 
1. The last case corresponds to electron- ion equiparti- 
tion. We also calculate spectra both for non-equilibrium 
ionization (NEI) models and for collisional ionization 
equilibrium (CIE) models for comparison. For the NEI 
models, the results for the model with the small shock 
electron heating efficiency (3 — 1/1800 will be discussed 
extensively throughout the paper, as this model maxi- 
mizes the departures from equilibrium in the outer re- 
gions of clusters. For the CIE models, we present a non- 
equipartition model with 13 ~ 1/1800 (CIE-Non-Eq) and 
an equipartition model with (3=1 (CIE-Eq). 

4.2. X-ray Spectra 

The projected rest frame spectra for several different 
models at two projected radii are shown in Figure [6l In 
the upper panels, we show the NEI model with a very 
small shock heating efficiency {/3 = 1/1800), and elec- 
trons and ions are in non-equipartition. The CIE-Non- 
Eq model is shown in the lower panels. The left panels 
show spectra at a radius of r = 2 Mpc, while the right 
panels show spectra at a radius of r = 3.5 Mpc. Each 
spectrum is binned with a bin size of Alog(i?) — 0.005. 

At r = 2 Mpc, the overall spectra are dominated by 
the free-free continuum emission over a wide range of en- 
ergy for both the NEI and the CIE-Non-Eq models. The 
continuum spectra for both models are nearly identical 
because of the dominant free- free emission with the same 
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electron temperature. The line emission is also very sim- 
ilar for both models because at this radius, the ionization 
timescale parameter is rather large (r ~ 8 x 10^^ cm~^ s). 
The most notable difference is the line intensity of the 
O VII triplet lines near ^ 0.57 eV. This line intensity for 
the NEI model is much higher than that of the CIE-Non- 
Eq model. There is almost no O VII line emission for the 
CIE-Non-Eq model. The strong O VII at r = 2 Mpc 
for the NEI model is mainly due to the projection of 
emission from the under-ionized O VII ions in the outer 
regions with much shorter ionization timescale parame- 
ters. The differences for other strong emission lines are 
much smaller between the two models. There is slightly 
more soft line emission below about 1 keV for the NEI 
model. 

At r = 3.5 Mpc, there are significant differences in the 
spectra between the NEI and the CIE-Non-Eq models. 
For the CIE-Non-Eq model, the continuum emission still 
dominates the overall spectrum; for the NEI model, the 
soft emission is dominated by lines. One of the most 
obvious differences in the line emission between the two 
models is again for the O VII triplet. The line intensity 
for the NEI model is much higher than that of the CIE- 
Non-Eq model. The O VII triplet is weak for the CIE- 
Non-Eq model. By inspecting a number of line ratios 
at different radii, we found that the ratio of the O VII 
and O VIII line intensities can be used as a diagnostic 
for the degree of ionization equilibrium, and this will be 
discussed in Section below. 

4.3. Surface Brightness Profiles 

The rest frame radial surface brightness profiles inte- 
grated over various energy bands for the outer regions 
of clusters are shown in Figure [T] The NEI models with 
(3 — 1/1800 are shown in thick lines in both the up- 
per left and upper right panels. The CIE-Non-Eq and 
the CIE-Eq models are shown as thin lines on the up- 
per left and upper right panels, respectively. The ra- 
tios of 5'NEi/5'ciE-Non-Eq and S'NEi/5'ciE-Eq are shown 
below the corresponding panels. Comparing the NEI 
and the CIE-Non-Eq models tells us the effects of non- 
equilibrium ionization alone, while comparing the NEI 
and the CIE-Eq model tells us the total effects of both 
non-equilibrium ionization and non-equipartition. 

From the left panels of Figure [3 we can see that 
non-equilibrium ionization significantly enhances the soft 
(0.3-1.0 kcV) emission in the outer regions. For the 
NEI model, the soft emission has been increased by more 
than 20% at around 3 Mpc, and up to nearly an order 
of magnitude around the shock radius compared to the 
CIE-Non-Eq model. The increase of the soft emission is 
due to the line emission by the under-ionized ions. For 
the CIE-Non-Eq model, the surface brightness profiles 
in all energy bands shown in Figure [7] decrease rapidly 
out to the shock radius. For the NEI model, the de- 
crease in surface brightness in the soft band as a func- 
tion of radius slows down near ^ 3 Mpc, and the soft 
band surface brightness actually increases with radius 
from ^3.7 Mpc out to nearly the shock radius, where the 
surface brightness drops rapidly. Within about 2.3 Mpc, 
non-equilibrium ionization effect is less than 5% in the 
soft band. The non-equilibrium ionization effect on the 
medium (1.0-2.0 keV) band is not as dramatic as the 
soft band, and the maximum increase is about 70% near 



3.8 Mpc. The non-equilibrium ionization effect on the 
hard (2.0-10.0 keV) band is less than 5% in most re- 
gions, and the effect is to lower the hard emission near 
the shock radius. The soft emission dominates the overall 
X-ray band (0.3-10.0 keV) outside of ^3 Mpc, and the 
overall X-ray emission decreases more slowly than for the 
CIE-Non-Eq model beyond that radius. Near the shock 
radius, the surface brightness in the overall X-ray band 
for the NEI model is about a fact of 6 higher than that 
for the CIE-Non-Eq model. 

The right panels of Figure [7] show that in addition to 
the non-equilibrium ionization effect, non-equipartition 
will increase the soft emission by a significant factor near 
the shock radius. This occurs because in the CIE-Non- 
Eq and NEI models, the electron temperatures in the 
outer regions are lower than for the CIE-Eq model, and 
this also leads to more soft X-ray line emission. The 
surface brightness profile in the overall X-ray band for the 
NEI model is a factor of 5 higher than that for the CIE- 
Eq model near the shock radius. A detailed discussion 
on the difference betwe en the CIE-Non-Eq and CIE-Eq 
models can be found in lWong fc SarazinI pOOl ). 

4.4. O VII and O VIII Line Ratio 

The most prominent non-equilibrium ionization signa- 
ture in the X-ray lines are for the line ratio of O VII 
and O VIII. Figure [8] shows the spectra for models at 
r ~ 2 Mpc. The spectra are shown in the 0.5-0.7 keV 
range which covers the rest frame energies of the O VII 
and O VIII lines. The upper panel shows the NEI model 
with a very low electron heating efficiency /3 — 1/1800. 
The middle panel shows the NEI model with an interme- 
diate P = 0.5. The lower panel shows the CIE-Non-Eq 
model with (3 — 1/1800. Each spectrum is binned with 
a bin size of AE ^ 0.0001 keV. 

In the upper panel of Figure 13 the most prominent 
lines are the He-like O VII triplets at 561.0, 569.6, 
and 574.0 eV, the H-like O VIII doublets at 653.5 and 
653.7 keV, and the He-like O VII line at 665.6 eV. Here, 
we focus on the line ratio between the He-like O VII 
triplets and the H-like O VIII doublets as a diagnostic for 
non-equilibrium ionization. These two lines have been 
used to search for the WHIM as well as to study the 
non-equilibrium ionization of th e WHIM ()Cen fc Fand 
l2006t lYoshikawa fc Sasakil [20061) . We find that the He- 
like O VII and the H-like O VIII lines also show strong 
signatures of non-equilibrium ionization in the outer re- 
gions of clusters. 

The upper and lower panels of Figure [S] show that 
non-equilibrium ionization strongly enhances the He-like 
O VII triplets compared to the CIE-Non-Eq model. The 
H-like O VIII doublets for the NEI model are only slightly 
stronger than for the CIE-Non-Eq model. This suggests 
that the ratio of the O VII and O VIII lines can be used 
as a diagnostic for non-equilibrium ionization. The O VII 
and O VIII lines are similarly strong for the NEI mod- 
els with electron heating efficiencies (3 — 0.5 and 1/1800, 
and we suggest that the O VII and O VIII line ratio is a 
good diagnostic for the non-equilibrium ionization for a 
wide range of electron heating efficiencies. 

Figure [S] shows the surface brightness for the O VII and 
O VIII lines for the NEI model with ^ = 1/1800. The sur- 
face brightness of the lines are calculated by subtracting 
the continuum surface brightness from the total surface 
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brightness within narrow energy ranges of 556.0-579.0 eV 
and 648.5-658.7 eV for the O VII triplets and the O VIII 
doublets, respectively. The energy bands were chosen 
to cover the O VII and the O VIII lines with a spectral 
resolution of 10 eV which is the expected value for the 
outer arrays of the X-ray Microcalorimeter Spectrometer 
(XMS) on the IXO. The continuum surface brightness 
were calculated by fitting the spectrum to a power law 
model in the energy range of 0.5-0.7 keV, excluding the 
lines. This simulates the techniques likely to be used to 
analyze real observations. We also show the continuum 
surface brightness within the narrow energy bands used 
to extract the line surface brightness in Figure [HI 

The surface brightness of the O VII triplets is nearly 
constant from ^ 1 Mpc to ~ 3 Mpc, and then rises grad- 
ually. Note that a flat surface brightness at inner radii 
and rising surface brightness at larger radii is the signa- 
ture of a shell of emission at large radii seen in projec- 
tion. That is, most of the O VII emission is actually at 
large radii where the ionization time scale is short. Be- 
yond about 3.9 Mpc, the O VII surface brightness rises 
rapidly to a peak value, and then drops. The continuum 
in the 556.0-579.0 eV energy band drops from a radius 
of 1 Mpc out to the shock radius. Beyond ~ 2.8 Mpc, 
the O VII line emission dominates over the continuum in 
the 556.0-579.0 eV energy band. For the O VIII doublets, 
the surface brightness drops from a radius of 1 Mpc out to 
~ 3 Mpc, and then rises to a peak value at 4 Mpc. The 
O VIII surface brightness then drops beyond ~ 4 Mpc. 
The O VIII line emission dominates over the continuum 
in the 648.5-658.7 eV energy band beyond ~ 2 Mpc. 

Figure [TU] shows the line ratios of O VII and O VIII, 
S{0 VIII)/S'(0 VII), for different models. To com- 
pare the effect of non-equilibrium ionization alone, we 
can compare the line ratios between the NEI model 
(solid line) and the CIE-Non-Eq model (dashed line), 
while both models assume an electron heating efficiency 
(3 = 1/1800. Both models use the same non-equipartition 
electron temperature to calculate the spectra, but one of 
them assumes equilibrium ionization while the other one 
does not. At ^ 1 Mpc, the line ratio for the NEI model 
is more than a factor of 2 lower than that of the CIE- 
Non-Eq model. The difference increases as the radius 
increases, and the differences are over an order of magni- 
tude for radii beyond ~ 2 Mpc. For the CIE-Eq model 
(dotted line) , the line ratio is very similar to that of the 
CIE-Non-Eq model, except near the shock regions where 
the line ratio for the the CIE-Eq model is a factor of a 
few higher than the CIE-Non-Eq model. Both the line 
ratios of the CIE-Eq and the CIE-Non-Eq models are 
above 10 in most regions between 1-4 Mpc. NEI models 
with electron heating efficiencies /3 = 0.5 and 1.0 are also 
shown in Figure [TOl The effect of increasing the electron 
heating efficiency is to raise the electron temperature at 
the shock, and hence increase the ionization rates. From 
Figurefini we can see that increasing the electron heating 
efficiency only affects the line ratio by less than a factor 
of two in most regions shown compared to the NEI model 
with 13 = 1/1800. The line ratios for all the NEI models 
with different electron heating efficiencies are less than 
10 for radii beyond about 1.3 Mpc. In summary, the line 
ratios in the outer regions for all the NEI models we cal- 
culated are significantly smaller than those for the CIE 
models, and the differences are larger than an order of 



magnitude for most regions beyond '-^ 2 Mpc. Such large 
differences can be used to distinguish between NEI and 
CEI models in real observations. 

5. DETECTABILITY OF O VII AND O VIII LINES WITH 
IXO AND TESTING THE NON-EQUILIBRIUM 
IONIZATION EFFECT 

In this section, we estimate whether the O VII and 
O VIII lines in cluster outer regions and the non- 
equilibrium ionization signatures can be detected with 
IXO. We consider the NEI model with /3 = 1/1800 for a 
cluster with an accreted mass of Mgh = 1.53 x 10^^ Mq 
at low redshift. The choice is justified by the fact that 
the non-equilibrium ionization effect does not depend 
strongly on mass and the effect is larger at lower redshift 
(Section [3|), and that the surface brightness for massive 
and low redshift clusters is higher. 

The X-ray Microcalorimeter Spectrometer (XMS) 
planned for the IXO can potentially detect the O VII 
and O VIII lines in cluster outer regions. The Wide Field 
Imager (WFI) does not have enough spectral resolution 
(> 50 eV) at around 0.6 keV, and the O VII triplets 
cannot be resolved from the 0.5 keV nitrogen line. The 
X-ray Grating Spectrometer covers the interesting en- 
ergy range, but the collecting area is too small to detect 
the weak O VII and O VIII lines. Therefore, we only con- 
sider the XMS in our estimations. The XMS has inner 
(core) and outer microcalorimeter arrays with expected 
spectral resolutions of 2.5 and 10 eV, respectively. We 
consider two cases when observing with the XMS. The 
first case (XMSC) is that only the inner core array is 
used, and the second case (XMSF) is that the full ar- 
ray (both inner and outer arrays) is used. For simplicity, 
when considering the XMSF, we assume the spectral res- 
olution to be the same as the outer array throughout the 
full array. This doesn't strongly affect the detectability 
of the lines, since the bands used to determine the fluxes 
in the lines are set by the line width in the outer array, 
and the line width in the inner core doesn't affect the line 
flux. The expected effective areas {A^s) for the XMS is 
about 10000 cm^ at around 0.6 keV. The relevant in- 
strument parameters for the two cases we considered are 
fisted in Table [1] 

5.1. Backgrounds 

The major background for the IXO observations are 
Non- X-ray Background (NXB), the soft emission from 
the local Galactic background (GXB), and these are in- 
cluded in our signal-to-noise ratio calculations. The clus- 
ter continuum emission will be much weaker than the line 
emission in the band widths we are interested in, but we 
also include the cluster continuum emission in our calcu- 
lations. With the spatial resolution of 5" and the very 
long exposure time needed for the line detection, point 
source contaminations will be negligible, and hence it is 
not included in our calculations. In fact, the GXB we 
used has included a component from unresolved AGNs, 
and this may overestimate the total background. The 
total count rates of the backgrounds we used for the two 
instrument setups are listed in Table [5J and are discussed 
below. 

We use the GXB simulated bv lFang et al] (|2005[ ) which 
included two thermal components to represent the Lo- 
cal Hot Bubble emission and the transabsorption emis- 
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sion ()SnowdeD|[l99l iKuntz fc Snowdenl I2000D . and one 
continuum emission component to represent unresolved 
AGN background. The para meters used for the i r back - 
ground model are based on iMcCammon et aTl (|2002| ). 
The most prominent emission around 0.6 keV is the 
lin e emission from nitrogen and oxygen ions (Figure 6 
in iFang et aIll2005D . With the very high XMS spec- 
tral resolutions, the O VII and O VIII lines from clus- 
ters beyond z k. 0.028 should be separated from the 
strong line emission in the GXB (see Table [3] below). 
We adopt the continuum intensity value of /gxb — 24 
photons cr n~^ s~^ sr~ ^ keV" -^ at around 0.6 keV, which 
is used in iFang et aTl ()2005l ) . The total count rate for 
the whole field of view (ripov) which covers the energy 
range of the lines is given by -Rgxb = -^gxb x ^cff x 
f^FOV X Aii'band, whcrc AiJband is the bandwidth cov- 
ered the lines. For the O VIII doublets, since the line 
separation is much smaller than the spectral resolutions 
of all the instruments we considered, AE^band is simply 
the spectral resolution of the corresponding instrument. 
For the O VII triplets, the three lines should be well sep- 
arated by the XMS core array. Hence, AE'band of O VII 
is then three times the spectral resolution of the XMS 
core array. However, for the XMS full array, the O VII 
triplets cannot be resolved. Therefore, Ai?band is then 
given by the maximum separation of the triplets (13 eV) 
plus the corresponding spectral resolutions. 

For a future mission like the 1X0^ the NXB is rather 
uncertain. We use the count rate of Fnxb = 8.1 x 10^^ 
photons s~^ arcmin"^ keV~^ at 0.6 keV for th e XMS 
estimated by the IXO team ()Smith et all I2O1O0 . The 
total count rate for the whole field of view which cov- 
ers the energy range of the lines is given by i?NXB = 
-Fnxb X i^FOV X A_Eband- To address the effects of the 
uncertainties, we also multiply the NXB by factors of 0.5 
and 2 in our calculations (Section 15.21 below) . 

The cluster continuum emission within the very nar- 
row energy bands we are interested in are much weaker 
than the line emission in cluster outer regions (Figure [9]) 
and the GXB. To detect the O VII and O VIII lines, it 
is also important to observe regions where the cluster 
continuum emission is weak compared to the line emis- 
sion. Therefore, the cluster continuum emission should 
not be important when estimating the signal-to-noise 
ratio. Nevertheless, we have included the cluster con- 
tinuum emission in our calculations. To be conserva- 
tive, we simply take the cluster continuum emission at 
r = 2.8 Mpc where the cluster continuum emission across 
648.5-658.7 eV (energy range where the O VII triplets 
are covered by the XMS full array) equal to the O VII 
emission. The continuum count rates -Rcont are listed in 
Table El 

5.2. Signal-to-Noise Ratio 

In order for the lines to be separated from the local 
GXB, clusters to be observed should be at high enough 
redshifts. Since the O VII triplets spread across 13 eV, 
the best targets should have redshifts such that the lines 
can be shifted by at least 13 eV plus the spectral resolu- 
tion. The targets should also not to have too high red- 
shift such that the O VII triplets will not overlap with 
the 0.5 keV nitrogen line. Table [3] lists the minimum 
and maximum redshifts for clusters to be observed by 
the XMS core and the XMS full arrays which meet these 



requirements. 

To calculate the signal of the O VII and O VIII lines for 
clusters for both the XMSC and XMCF cases, we assume 
the nominal massive cluster model to be at a redshift of 
0.05. The redshift is chosen such that it is slightly higher 
than Zmin for the XMS full array in Table El It is also 
important that there actually be clusters which are at 
or within the selected redshift which are fairly regular 
in shape in their outer regions and with masses which 
are comparable to the cluster model we calculated (e.g., 
Abell 85, Abeh 1795). 

To calculate the total count rates of the lines from 
a cluster at the assumed redshift as observed by IXO 
(^line), we first convolve the rest frame surface bright- 
ness of the lines (Figure [9]) with a top hat function 
with a width of 0.114 (0.286) Mpc which corresponds 
to the physical distance at z = 0.05 covered by the 
XMSC (XMSF) field of view, and then muhiply the 
convolved rest frame surface brightness by a factor of 
AeffriFOv(47r)~^(l -I- z)~^ to give the count rate i?iino- 
In calculating the signal-to-noise ratios, count rates at a 
radius where the surface brightness of the weaker line is 
maximum in the outer regions are used as the optimum 
model count rates. The optimum count rates are also 
listed in Table d 

We calculate the signal-to-noise ratios for the optimum 
models (1 NXB, 1 GXB) which use the count rates listed 
in Table [21 The signal for each line is given by ^ctn = 
^linc t, where t is the exposure time. The noise for each 
line is given by 

N — \J (i?iinc + ^cont + -Rgxb + ^nxb) t ■ (18) 
The signal-to-noise ratio of the line ratio is given by 
SN{0 VIH/0 VH) = 

[[SN{0 VIII))-2 + (^SN{0 VII))-2]-i/2 ^ ^^g^ 

where SN{0 VII) and SN{0 VIII) are the signal-to- 
noise ratios of the O VII and O VIII lines, respectively. 

To address the effects of the NXB uncertainties, we 
also calculate the signal-to-noise ratios for models with 
the NXB muhiplied by a factor of 2 (2 NXB, 1 GXB) and 
a factor of 0.5 (0.5 NXB, 1 GXB). The GXB varies with 
sky position. To address the uncertainties in the GXB, 
we vary the GXB by a factor of 0.8 (1 NXB, 0.8 GXB), 
1.2 (1 NXB, 1.2 GXB), and 1.5 (1 NXB, 1.5 GXB). In 
real observations, the line signals may not be optimum. 
We consider models with the count rates of the lines to 
be a factor of 0.75 (0.75 i?iine), 0.50 (0.50 i^uno), and 0.25 
(0.25 -Riinc) of the optimum models. 

The signal-to-noise ratios as a function of exposure 
time for the different models are shown in Figure 1111 
To ensure that there are enough counts for the lines, we 
only plot the signal-to-noise ratios if the total count for 
each line S'ctn is larger than 30. 

For the XMSC instrument setup with the optimum 
model (1 NXB, 1 GXB), in order to get a signal-to-noise 
ratio of 3 for the O VII triplets, about 220 ksec is needed. 
With a deeper observation of 600 ksec, the signal-to-noise 
ratio can exceed 5. For the O VIII doublets, a shorter ex- 
posure time of ^ 180 ksec is need to have S'ctn > 30 with 
a signal-to-noise ratio of 3.2. About 450 ksec is needed 
to get a signal-to- noise ratio of 5, and it is possible to get 
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a signal-to-noise ratio up to 7 with a very deep observa- 
tion of 900 ksec. If no single cluster is observed for such 
a long exposure, the very deep exposure can be achieved 
by stacking many observations of outer regions of many 
clusters. In this case, the line emission measured is av- 
erage over many clusters. For a given exposure time, in- 
creasing the NXB by a factor of 2 or increasing the GXB 
by a factor of 1.5 only decrease the signal-to-noise ratios 
by about 10% for both the O VII or O VIII hues. As ex- 
pected, decreasing the line signals significantly decreases 
the signal-to-noisc ratios. If the line signals are smaller 
than half of the optimum values, it will be impossible to 
detect the O VII (O VIII) line with an exposure time less 
than - 0.8 (0.5) Msec. Limited by the O VII detection, 
an exposure time of about 220 ksec is sufficient to have 
a 2. 3- (7 determination for the line ratio for our optimum 
model. A 3-a measurement of the line ratio will require 
an exposure time of about 380 ksec. 

If the electron temperature can be measured from ob- 
servations such as X-ray spectroscopy or hardness ratios, 
the measured line ratio can be compared to the CIE ratio 
inferred by the electron temperature, and this will pro- 
vide a direct test to the ionization states of the plasma 
without ambiguity. In the outermost regions of a cluster, 
the surface brightest may be too low and hence reliable 
temperature measurement may not be feasible. However, 
as shown in Figure [TOl beyond a radius of ^ 2.5 Mpc 
where the electron temperatures are always higher than 
about 0.5 keV for either the non-equipartition model 
with f3 = 1/1800 or the equipartition model, the CIE 
line ratios are always higher than ^4. In contrast, the 
NEI line ratios are always lower than ^ 2 for the elec- 
tron temperature range of interest in realistic clusters 
for all the models we considered. Hence, a 3-ct measure- 
ment of the line ratio < 2 (A[S'(0 VIII)/5(0 VII)] < 
2/3 at 1-a) is a strong evidence that the ions are in 
non-equilibrium ionization due to the huge difference be- 
tween the NEI and CIE line ratios for the temperature 
range of realistic clusters. If the line ratio is measured 
to be even lower (e.g., < 1), only a 2 or 3-(t will be suf- 
ficient to rule out the CIE model. On the other hand, 
if the line ratio is measured to be as high as 4 in the 
outermost regions, a 6-ct (A[S'(0 VIII)/5'(0 VII)] = 2/3 
at 1-a) will be necessary to rule out the NEI low line 
ratio of '--^ 2 at a 3-a level. Even if the line ratio is 
measured within about < 2 Mpc where the line ratio 
for the NEI model can be as high as ~ 4, the surface 
brightness at this relative small radius is high enough so 
that electron temper ature can b e measured by the IXO 
or even Suzaku (e.g.. iGeorge et a l. 2009; Rci prich et al.l 
[200l iBasu et al.l [20101 iHoshino'et al.i 1201 0) with confi- 
dence. To study whether the IXO XMS can measure 
the electron temperature at this radius, we have carried 
out simulations using XSPEC with our model spectrum 
taken at radius ^ 2 Mpc. The detector response and 
background files were taken from the latest IXO simu- 
lator SIMX0. Both the NXB and the GXB have been 
included in the background file. We assumed an absorp- 
tion model wahs with nn = 5 x 10^° cm~^. We found 
that the IXO XMS core detector can collect about 2400 
photon counts (0.2-3.0 keV) for a 150 ksec observation. 
When fitting the continuum of the spectrum, we assumed 

^ http:/ /ixo. gsfc.nasa.gov/scicncc/simulator. html 



an absorbed thermal model {wabs*apec) and ignored the 
strong lines. Using the wabs*nei model, the temperature 
determined is within the uncertainty of the wabs*apec 
model and does not affect the results significantly. We 
found that the electron temperature can be constrained 
to Te = 4.0t} j keV. The CIE line ratios within the de- 
termined temperature range are always higher than 13, 
and hence a 2 to 3-cr measurement for the low line ratio 
will also be enough to test the non-equilibrium ioniza- 
tion model. In fact, the IXO WFI should constrain the 
electron tem perature better be cause of its lower particle 
background ([Smith et al.ll2010[) . 

For the XMSF instrument setup, with the optimum 
model (1 NXB, 1 GXB), about - 130 (100) ksec is needed 
to get a signal-to-noise ratio of 3 for the O VII triplets 
(O VIII doublets). A longer exposure of ^ 350 (270) ksec 
is needed to get a signal-to-noise ratio of 5 for the O VII 
(O VIII) lines, and a very deep exposure of ~ 1 (1) Msec 
is needed for a 8 (lO)-cr detection. Increasing the NXB 
by a factor of 2 or increasing the GXB by a factor of 
1.5 decreases the signal-to-noise ratios by less than ~ 
13% for both the O VII and O VIII hues. Similarly to 
the XMSC instrument setup, the signal-to-noise ratios 
decrease as the line signals decrease. For the optimal 
model, a shorter exposure time of about 130 ksec (also 
limited by the O VII detection) compared to the XMSC 
setup will give the same 2. 3-ct detection for the line ratio. 
An exposure time of about 230 ksec will be needed for a 
3-CT measurement of the line ratio. 

Overall, the XMSF instrument setup can achieve 
higher signal-to-noise ratios for a given exposure time 
compared to the XMSC setup for the optimum mod- 
els. This is mainly due to the larger field of view of 
the detector which can collect more photons with the 
same exposure time. The XMSF is slightly more sub- 
ject to the continuum background uncertainties (cluster 
continuum emission, NXB, and GXB) due to the poorer 
spectral resolution. For both the XMSC and XMSF se- 
tups, detecting the O VII and O VIII lines and testing 
non-equilibrium ionization in cluster outer regions are 
promising. 

6. DISCUSSIONS AND CONCLUSIONS 

Studying the physics in the outer regions of clusters is 
very important to understand how clusters are formed, 
how the the intracluster gas is heated, as well as to con- 
strain the formation of large-scale structure. Because of 
the very low density in cluster outer regions, the col- 
lisional timescales are very long and comparable to the 
cluster age. Electrons and ions passed through the accre- 
tion shocks may not have enough time to reach equipar- 
tition and the ions may be un der-ionized. In a previous 
paper fW ong fc Saraa[200g ). we have studied the non- 
equipartition effects on clusters using one-dimensional 
hydrodynamic simulations. In this paper, we systemat- 
ically studied non-equilibrium ionization effects on clus- 
ters and the X-rays signatur es using the same set o f sim- 
ulations we have developed (iWong fc Sara zin 2009) . 

By using semi-analytic arguments together with nu- 
merical simulations, we have shown that the non- 
equilibrium ionization effect is nearly independent of 
cluster mass but depends strongly on redshift. In partic- 
ular, non-equilibrium ionization effects are stronger for 
low-redshift clusters. Therefore, the brighter massive 
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clusters at low-rcdshifts are good candidates for studying 
the non-equilibrium ionization effects. 

We systematically studied non-equilibrium ionization 
signatures in X-rays for a massive cluster with Mgh = 
1.53 X 10^^ Mq. We first calculated the ionization frac- 
tions for 11 elements heavier than He following the elec- 
tron temperature and density evohitions of each fluid el- 
ement. We then calculated the X-ray emissivity of each 
fluid element and the resulting projected spectra for the 
cluster. Since the electron temperature profiles depend 
on electron heating efficiency /3, we have considered three 
different possibilities which represent a very low heat- 
ing efficiency {(3 = 1/1800), an intermediate heating ef- 
ficiency {f3 — 0.5), and equipartition, /? = 1. We also 
considered models which assume equilibrium ionization 
for comparison. 

At a radius (e.g., 2 Mpc) where the ionization timescale 
is long, the overall spectra for the NEI and CIE-Non-Eq 
models are very similar. This is because of the dom- 
inant free-free emission, and both models assume the 
same electron temperature. However, in the outer re- 
gions, e.g, at r ~ 3.5 Mpc which is between iJyir and 
Rsh, the soft emission in the NEI model is dominated by 
line emission, where the CIE-Non-Eq spectrum is still 
dominated by the contimuim free-free emission. 

By analyzing the surface brightness profiles, we found 
that soft emission (0.3-1.0 keV) for the NEI model can be 
enhanced by more than 20% at around 3 Mpc, and up to 
nearly an order of magnitude near the shock radius com- 
pared to the CIE-Non-Eq model. The soft emission en- 
hancement is mainly diia to the line emission from under- 
ionized ions. The non-equilibrium ionization effects on 
the medium (1.0-2.0 keV) and hard (2.0-10.0 keV) band 
emissions are smaller. The overall X-ray band (0.3- 
10.0 keV) emission is dominated by the soft emission, 
and the total X-ray emission for the NEI model decrease 
much slower than that of the CIE-Non-Eq model. Thus, 
if cluster outer regions are in non-equilibrium ionization, 
the shock region will be much more luminous compared 
to the CIE-Non-Eq model. 

By inspecting a number of spectra, we found that the 
most prominent non-equilibrium ionization signature in 
line emission is the line ratio of the He-like O VII triplets 
and the H-like O VIII doublets, 5'(0 VIII)/S'(0 VII). 
The line ratios for the CIE models are higher than 10 for 
most regions between r = 1-4 Mpc, while the line ratios 
are smaller than 10 for the NEI models. The differences 
in the line ratios between the NEI and CIE models in- 
crease with radius, and the differences are more than an 
order of magnitude for radii beyond ~ 2 Mpc. These 
results are insensitive to the degree of non-equipartition 
or electron heating efficiency /?. We suggest that the line 
ratios can be used to distinguish between the NEI and 
CIE models. The electron temperature profile can be de- 
termined from fits to the continuum spectra of the outer 
regions of clusters, allowing the CIE line ratios to be 
determined. Comparison to the observed ratios should 
show the effects of non-equilibrium ionization. Note that 
a line ratio of 5(0 VIII) /S{0 VII) < 3 in the outer region 
of a massive clusters is a clear signal of NEI. 

We have also studied the detectability of the O VII 
and O VIII lines around cluster accretion shock regions 
with IXO, as well as the test for non-equilibrium ioniza- 
tion using the line ratio. For our optimum model, we 



found that with the XMS core array, an exposure time 
of 220 ksec is need to have a S.O-ct detection of the O VII 
lines and about 180 ksec is need to have > 30 counts for 
a 3.2-cr detection of the O VIII lines. The uncertainties 
in NXB and GXB will not afi^ect the results significantly. 
For the XMS full array while we assume the spectral res- 
olution to be the same as the outer array throughout the 
detector, we found that the signal-to-noise ratios for our 
optimum model are higher for the same exposure time 
as the XMS core array. In particular, only about 130 
(100) ksec is needed to detect the O VII (O VIII) line. 
The XMS full array is only slightly more subject to NXB 
and GXB imcertainties due to the poorer spectral reso- 
lution. 

To test the non-equilibrium ionization model without 
ambiguity requires measurements of both the electron 
temperature (or hardness ratio) and the line ratio so 
that the measured line ratio can be compared to the CIE 
line ratio inferred by the electron temperature. We have 
shown that this can be done within < 2 Mpc of a cluster 
by the IXO with sufficient confidence. Beyond ~ 2.5 Mpc 
where the surface brightness may be too low and mea- 
suring the electron temperature may be diSicult, we have 
shown that if the line ratio is measured to be as low as 
~ 2 at 3-cr (A[S'(0 VIII)/S'(0 VII)] - 2/3 at 1-cr) in the 
outermost regions, this will rule out the CIE model at a 
3-a level since the CIE line ratios are always higher than 
4 for realistic cluster temperatures. A 3 or 4-c7 measure- 
ment of such a low line ratio is sufficient to provide a 
strong test of the non-equilibrium ionization. If the line 
ratio is measured to be even lower (e.g., < 1), only a 2 
or 3-(7 will be sufficient to rule out the CIE model. On 
the other hand, if the line ratio is measured to be as 
high as 4, a 6-cr measurement (A[S'(0 VIII)/S'(0 VII)] 
= 2/3 at l-cr) will be necessary to rule out the NEI low 
line ratio of ^ 2 at a 3-(t level. We found that an ob- 
servation with about 130 (220) ksec with the XMS full 
(core) array is enough to measure the line ratio at 2.3- 
fj. For a 3-cr measurement of the line ratio, about 230 
(380) ksec will be needed for the XMS full (core) array, 
and this will provide a strong test for non-equilibrium 
ionization. In summary, detecting the O VII and O VIII 
lines around the cluster accretion shock regions and test- 
ing non-equilibrium ionization in cluster outer regions 
with IXO are promising. 

It is expected that the O VII and O VIII lines from 
WHIM will also be strong. Because of the high spec- 
tral resolution of the XMS, emissions from different red- 
shifts should be easily separated. Only the emission 
from WHIM immediately surrounding the target clus- 
ter will be potentially confused with the emission from 
cluster outer regions. To observe the O VII and O VIII 
lines and study NEI effects in cluster accretion shock re- 
gions, it will be best to avoid observing directions along 
the filaments where it is believed that denser preheated 
WHIMs and subclusters are preferentially accreted onto 
more massive and relaxed galaxy clusters. 

What do we learn about clusters from the ioniza- 
tion state of the outer gas? Since collisional ionization 
and recombination rates involve straightforward atomic 
physics, the processes are not in question and the rates 
are reasonably well-known. Unlike shock electron heat- 
ing or rates for transport processes like thermal conduc- 
tion, the basic physics is not uncertain and magnetic 
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fields do not affect the results in a significant way. What 
we mainly learn about is the pre-shock physical state of 
materials which are being accreted by the cluster. If most 
of the WHIM is ionized beyond O VII, then the effects 
described in this paper will be greatly reduced. If most of 
the material currently being accreted by clusters comes 
in through filaments which have a higher ionization, then 
NEI effects will be diminished significantly. If most of the 
gas being added to clusters at present comes in through 
mergers with groups which deposit most of the gas in the 
inner regions of clusters, the gas will achieve CIE quickly. 

From the theoretical point of view, with the increasing 
number of observations of galaxy cluster outer regions 
(~ i?20o) and the potential to extend observations out 
to the shock radius with IXO in the future, it is neces- 
sary to perform more detailed simulations than ours. It 
is also interesting to extend our work to study the con- 
nections between the shocked ICM and the more diffuse 
WHIM surrounding clusters. Three-dimensional simula- 
tions will be essential to understand the effects of merg- 
ers or filament accretion on the degree of ionizations in 
different regions of clusters. This will allow us to char- 
acterize the variation of non-equilibrium signatures in 
the clusters; such calculations are essential to compare 
observational signatures with our understanding of the 
cluster physics near the accretion shocks. Cosmologi- 
cal simulations have been performed recentl y to study 
the NEI signatures (,Yoshikawa et al.,.2003.: Cen fc Fang! 



I2006t lYoshikawa fc Sasakil I2OO60 . These studies have 
shown that both non-equipartition and non-equilibrium 
ionization effects are important in cluster outer regions; 
although they focus more on the lower density and lower 
temperature WHIM. High resolution simulations were 
also performed for studying NEI effects in clusters, but 
these are limited to binary mergers with idealized ini- 
tial conditions and focus on the denser merger shocks 
(|Akahori &: Yoshikawal[2008ll2010[ ). Re-simulating repre- 
sentative clusters and the surrounding WHIM from cos- 
mological simulations with higher resolutions and includ- 
ing realistic physics (e.g., cooling, conduction, turbulent 
pressure, magnetic pressure, and relativistic support by 
cosmic rays) will be necessary to provide realistic model 
images and spectra. The different observational signa- 
tures and connections between the ICM and the more 
diffuse WHIM can also be addressed self-consistently by 
these simulations. 
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TABLE 1 

Parameters for the IXO XMS Core and Full 
Arrays 



Instrument 


Aeff (cm2) 


fipov (arcmin^) 


AE (eV) 


XMSC 


10000 


2x2 


2.5 


XRISF 


10000 


5 X o 


10 



Note. — Column 3 lists the field of view. Column 
4 lists the spectral resolution (FWHM). 



TABLE 2 
Count Rates 



Instrument 


Ion 


-^line 


-Rcont 










(10-* counts s-i) 


(10-" counts s^i) 


(10-4 counts s-i) 


(10-* counts s-i) 


XMSC 


o VII 


2.11 


0.210 


2.44 


6.09 


XMSF 


O VII 


11.3 


4.03 


46.7 


117 


XMSC 


o VIII 


1.66 


0.171 


0.878 


2.19 


XMSF 


o VIII 


8.69 


1.43 


20.8 


51.8 



TABLE 3 

MiNiirrM ANn MAXiirrM Redsiiifts 



Instrument 






^-£"111 ax 


■2^max 




(eV) 




(eV) 




XMSC 


15.5 


0.0278 


58.1 


0.116 


XMSF 


23.0 


0.0417 


50.6 


0.0992 



Note. — Column 2 lists the minimum en- 
ergy shifts in order for the O VII triplets to be 
separated from Galactic lines. Column 3 lists 
the corresponding rcdshifts of AiJmin- Column 
4 lists the maximum energy shifts in order for 
the O VII triplets not to overlap with the 0.5 keV 
nitrogen line. Column 5 lists the corresponding 
redshifts of ASmax- 
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Fig. 1. — Ionization timescale parameter r (Equation [2]l versus radius r scaled to the cluster shock radius -Rgh for cluster models with 
total accreted masses of Msh = 0.77 (dashed line), 1.53 (sohd hne), and 3.06 (dotted line) X lO^^^ Mq at a redshift of 2 = 0. The shock 
radii for the three clusters from small to high mass are figji = 3.31,4.22, and 5.41 Mpc, respectively. All three lines lie almost on top of 
one another. 




Fig. 2. — Dimensionlcss ionization parameter Ci{Te)T for O VII versus the scaled radius {r/R^h) for clusters with total accreted masses 
of = 0.77 (dashed line), 1.53 (solid hne), and 3.06 (dotted line) X lO^^ jvfg at a redshift of z = 0. Models with /3 = 1/1800 

(non-equipartition) and /3 = 1 (equipartition) are shown in thick and thin lines, respectively. 
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Fig. 3. — Dimensionless ionization parameter Ci{T^) t for O VII versus the scaled radius {r/R^y^) for clusters at different redshifts. The 
cluster model with a total accreted mass M^h = 1-53 x 10^^ Mq at = is used. Models with /3 = 1/1800 (non-equipartition) and /3 = 1 
(equipartition) are shown in thick and thin lines, respectively. 




Fig. 4. — Non-equilibrium ionization bias versus the scaled radius (r/Rgh) for clusters with total accreted masses of Mgh = 0.77 (dashed 
Une), 1.53 (solid line), and 3.06 (dotted line) x 10^^ Mq at a redshift of 2 = 0. Models with /3 = 1/1800 (non-equipartition) and /3 = 1 
(equipartition) are shown in thick and thin lines, respectively. 
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Fig. 5. — Non-equilibrium ionization bias versus the scaled radius (r/R^^^) for clusters at different redsliifts. The cluster model with a 
total accreted mass Mgh = 1-53 x lO'^^ Mq at 2: = is used. Models with j3 = 1/1800 (non-equipartition) and /3 = 1 (equipartition) are 
shown in thick and thin lines, respectively. 



r = 2Mpc r = 3.5Mpc 




Energy (keV) 



Fig. 6. — Projected rest frame spectra for the cluster model with an accreted mass of Msh = 1.53 X 10^5 M0 at 2 = are shown at two 
projected radii (left: 2 Mpc; right: 3.5 Mpc). Upper panels: the NEI model with /3 = 1/1800 and non-equipartition of electrons and ions. 
Lower panels: projected rest frame spectra for the CIE model with the same /3 and non-equipartition. 
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r (Mpc) r (Mpc) 



Fig. 7. — Upper left panel: rest frame projected surface brightness profiles for different energy bands for the NEI model with /3 = 1/1800 
(thick lines) and the CIE-Non-Eq model (thin lines). Upper right panel: rest frame projected surface brightness profiles for the NEI 
model with /3 = 1/1800 (thick lines) and the CIE-Eq model (thin lines). Lower left panel: ratios of the surface brightness profiles 
•SNEl/SciE-Non-Eq- Lower right panel: ratios of the surface brightness profiles 5NEl/<SciE-Eq- 
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Fig. 8. — Upper panel: projected rest frame spectra for the NEI model with /3 = 1/1800. Middle panel: projected rest frame spectra for 
the NEI model with f} = 0.5. Lower panel: projected rest frame spectra for the CIE-Non-Eq model with /3 = 1/1800. All spectra are for 
the Msh = 1.53 X 10^5 M© at z = 0, and are for a projected radius of r = 2 Mpc. The spectra are binned with Afi = 0.1 eV. 
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Fig. 9. — Surface brightness of the O VII and O VIII hncs for the NEI model with /9 = 1/1800 as a function of projected radius. The 
continuum emission within the narrow energy bands covering the lines is also shown. 
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Fig. 10.— Line ratios 5(0 VIII)/S(0 VII) for the NEI models with /9 = 1/1800 (solid), 0.5 (dash-dotted), and 1.0 (dash-dot-dotted). 
Line ratios for the CIE-Non-Eq (dashed) and CIE-Eq (dotted) models are also shown. 
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Fig. 11. — Left (middle) panels show the signal-to-noise ratios for the O VII triplets (O VIII doublets) expected to be detected by IXO. 
The right panels show the signal-to-noise ratios for the O VII and O VIII line ratios. The upper panels correspond to the XMSC, and the 
lower panels correspond to the XMSF. Different signal and noise levels are assumed for the different models. The legends for the lines are 
the same for all panels, but are separated into two panels due to the space limitation. 



